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ABSTRACT 

Cosmic filaments play a role of bridges along which matter and gas accrete onto 
galaxies to trigger star formation and feed central black holes. Here we explore the 
correlations between the intrinsic properties of void galaxies and the linearity Rl oi 
void filaments (degree of filament's straightness) . We focus on void regions since the 
bridge effect of filaments should be most conspicuous in the pristine underdense regions 
like voids. Analyzing the Millennium-Run semi-analytic galaxy catalogue, we identify 
void filaments consisting of more than four galaxies (three edges) and calculate the 
means of central black hole mass, star formation rate, and stellar mass as a function of 
It is shown that the void galaxies constituting more straight filaments tend to have 
higher luminosity, more massive central black holes and higher star formation rate. 
Among the three properties, the central black hole mass is most strongly correlated 
with Rl. It is also shown that the dark halos constituting straight filaments tend to 
have similar masses. Our results suggest that the fuel-supply for central black holes 
and star formation of void galaxies occurs most efficiently along straight void filaments 
whose potential wells are generated by similar-mass dark halos. 
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1 INTRODUCTION 

The spatial distribution of galaxies in the universe is quite 
anisotropic, forming a filamentary web on large scales, which 
is often called the cosmic web. The existence of cosmic 
web has been explained as a natural phenomenon caused 
by the large-scale coherence in the primordial tidal field 
(|Bond et all Il996l '). The physical properties of the galax- 
ies located in the cosmic web are believed to be related 
to the merging/accretion history of the underlying dark 
halos, which in turn are affected not only by the local 
density but also by the tidal fields as manifested by the 
cosmic web. Although the effect of the tidal field on the 
galaxy's properties is expected to be weaker than that of 
the density field, several numerical and observational evi- 
dences have recently been reported for t he tidal influence s 
on dark halos and galaxies. For instance. iDesiacguesI (|2008| ) 
has shown by N-body simulations that the tidal effect from 
the surrounding matter distribution causes earlier virializa- 
tions of dark hal os embedded in overdense environments. 
iHahn et al.l l|2008l ) have claimed using the results from high- 
resolution N-body experiments that the high densities of 
the r egions with str ong tidal forces induce the assembly 
bias (|Gao et all 120051 ). The link between the tidal field and 
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galaxy properties have been also detected in real observa- 
tions. It has been recently shown that the galaxy morphol- 
ogy and luminosity depend on the l arge scale tidal field 
(|Lee fc Erdogdul2007l : lLee fc Leell2008l ). When the density is 
constrained to a narrow range, more luminous and late-type 
galaxies are found to reside in regions with higher ellipticity. 

In previous works, however, what has drawn little at- 
tention is the effect of the directions of the tidal forces on 
the galaxy properties. If the halo formation and the merg- 
ing/accretion events occur preferentially along filaments, 
i.e., in the directions of the principal axes of the tidal field, 
the tendency of the spatial alignments of the tidal field 
should be taken into account when the tidal infl uences on the 
galaxy properties are explored. Very recently, iPorter et al.l 
( 2008.) have reported that the star formation rate of a galaxy 
is observed to be enhanced when it falls along a superclus- 
ter fllament. This observational evidence suggests that the 
merging and accreting along the coherent directions of the 
tidal field trigger the star formation more efficiently. 

In practice, the degree of the tidal alignment can be 
quantified by measuring the straightness of a filament. The 
more straight a filament is, the stronger the tidal alignment 
is. Here, we propose a hypothesis that the highly coherent 
accretion of matter and flow of gas along straight filaments 
trigger the star formation and feed the central black holes 
most efficiently. We call it the bridge effect of the filaments. 
To test this hypothesis, we investigate the cross-correlations 
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between the stellar mass, central black hole mass, and star 
formation rate of the void galaxies and the straightness of 
the filaments where the galaxies are located. The reason that 
we consider only void regions for this investigation is that 
the strong effect of the density field can be controlled to a 
minimum level in the extreme underdense regions like voids. 

The organization of this paper is as follows. In §2, we 
analyze the data from the Millennium Run semi-analytic 
catalogue. In §3, we explore the correlations between the 
properties of void galaxies and the linearity of host filaments. 
A summary of the results and the final conclusion are given 
in 54. 



2 DATA 

We use the void catalogue that w as constructed in 
our previous work (IPark fc Led l2007all b y app lying the 
void-finding algorithm of lHovle' fc VogelevI (|2002l . hereafter 
HV02) to the Millennium - Run semi-analytic galaxy cata- 
loeu e (ISpring el ct al. 2005*'). As done in our companion pa- 
per l|Park fc Lee 2009. hereafter PL09), we first select those 
large voids which contain more than 30 galaxies and then 
apply the filament finding algorithm of I Barrow et al.1 ()l985l l 
to the select e d voi ds. The filament finding algorithm of 
iBarrow et al] (|l985l ') basically utilizes the minimal spanning 
tree (MST) technique to identify filamentary structures. A 
detailed description of the filament-identification in the Mil- 
lennium voids is provided in PL09. Here we provide a concise 
summary of it. 

For each selected void, we construct a MST of the void 
galaxies as follows: A starting galaxy (i.e., a node) is linked 
to its nearest neighbor galaxy by a straight line (i.e., an 
edge). This partial tree is extended by connecting it to the 
next nearest neighbor void galaxy. When all galaxies in a 
given void are connected, a MST is obtained. Then, we re- 
duce the MST of each void to identify the filaments through 
pruning and separating processes. Each MST is pruned at 
the p-level by removing the small-scale twigs. The pruned 
MSTs are separated into several distinct shorter pieces by 
removing those edges longer than a given length-threshold, 
Ic- The level of pruning and the length threshold are set at 
p = 5 and Ic — I + ai, respectively, where I and ai denote 
the mean edge-length averaged over all unreduced MSTs and 
its standard deviation, respectively (PL09). It is worth men- 
tioning that the pruning level was set at p = 4 in PL09 where 
the void filaments were found from the dark halos identified 
by the Friends-of- Friends (FOF) algorithm since the size dis- 
tribution of void filaments were found to become stabilized 
at p ^ 4. Now that the void filaments are found using not 
the FOF halos but the galaxies in the current analysis, we 
notice that the size distribution of void filaments becomes 
stabilized when the pruning level is p ^ 5. The origin of this 
difference in the value of p lies in the fact that the galax- 
ies correspond to the smaller-size substructures of the FOF 
halos. When the galaxies are used to construct the MSTs, 
the MSTs have more twigs and more superfluous small-scale 
noise. That is why a higher level p = 5 is required to prune 
the MSTs. 

The size of each filament, S, is measured as the spa- 
tial extent of the filament galaxies (see eq.[l] in PL09), and 
the linearity of a given filament Rl, which is introduced to 
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Figure 1. Mean sizes of the void filaments identified in the Mil- 
lennium Run semi-analytic galaxy catalogue as a function of the 
filament's linearity Rj^ for various edge-thresholds. 



quantify the filament straightness, is measure d as the ratio 
of th e end-to-end distance to the total length iBarrow et al.l 
1 198 5). The range oi Rl is [0, 1] and the degree of the filament 
straightness increases with Rl- Note that the value of Rl is 
biased toward the small value of S. The shortest filaments 
that contain only one edge (i.e., two nodes) have Rl = 1 by 
definition. To study the correlations between the properties 
of void galaxies and the linearity of void filaments, it should 
be required to select only long filaments since the short fila- 
ments are dominant in the high linearity section, causing a 
bias in the measurement of correlations. In other words, to 
find true correlations between Rl and void galaxy's prop- 
erties, it is necessary to put some lower limit {ric) on the 
number of filament's edges. Fig. [1] plots the mean sizes of 
the void filaments as a function of Rl, for the five different 
cases of the edge-threshold {uc = 1, 2, 3, 4 and 5 as dotted, 
dot-dashed, dashed, solid and long dashed line respectively). 
As can be seen, for all cases (5) decreases as Rl increases. 
But, for ric ^ 4 (long-dashed and solid lines), (S) does not 
decreases significantly with Rl in the high-linearity range 
of Rl ^0.7. Using this empirical result, we choose ric = 4 
and select only those filaments which consist of 4 edges or 
more (i.e., five galaxies or more). 

Table [T] lists the statistical properties of the void fila- 
ments. The first row corresponds to the case where no edge 
threshold is applied (one edge or more), while the second 
one corresponds to the case that only those filaments with 
He — 4 are considered. As can be seen, the mean linearity of 
the void filaments is low when the edge-threshold, ric = 4, is 
applied. It indicates that the short-size filaments dominate 
the high-linearity section. Fig. [2] illustrates two examples of 
the void filaments with two different values of Rl, demon- 
strating clearly that the higher the value of Rl is, the more 
straight a filament is. 



Bridging Effect of Void Filaments 3 




Figure 2. Illustration of the two examples of void filaments having two different values of the filament linearity, R^. 



Table 1. Edge-cut ric, total number of void filaments Nf, mean 
number of void filaments per a void (fij-), mean length of void 
filaments (Rf) in unit of h^^Mpc, and mean linearity of void 
filaments (Rl)- 
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3 THE BRIDGE EFFECT 

3.1 Linearity vs. Galaxy Property 

The void filaments correspond to deep potential wells within 
which excess of matter and gas can exist in otherwise ex- 
tremely underdense environments. They bridge the void 
galaxies onto which matter and gas from the surrounding 
large scale structures accrete and flow along their longest- 
axis directions. Therefore, the coherent accretion and flow of 
matter and gas may occur more efficiently in more straight 
filaments. It is expected that there may be a link between 
the physical properties of void galaxies and the linearity of 
their parent filaments, caused by this bridge effect of void 
filaments. 

We measure Rl of the selected void filaments and bin 
the values of their linearity Rl- Then, we calculate the 
mean values of central black hole mass (Mcbh), stellar mass 
(Afstar), total mass (Mtotai), and the star formation rate 
(SFR) of the void galaxies constituting a given void fila- 
ments with Rl in a, differential range of [Rl,Rl + cLRl]- 



Here MtotaEl means the sum of all gas masses in the cata- 
logue. Fig. 13] plots the means of the central black hole mass 
(top-left), stellar mass (top- right), total gas mass (bottom- 
left), and star formation rate (bottom-right) as a function 
of parent filament's linearity, Rl- In each panel the dotted 
line corresponds to the mean value averaged over all void fil- 
aments. The errors are calculated as (APq)^''^ /\^Nb where 
Nb is the number of void filaments belonging to a given bin 
of Rl and Pq represents a given galaxy property. As can 
be seen, the four different properties show a similar trend. 
The mean value of galaxy's property increases as the value 
of Rl increases. That is, the void galaxies constituting more 
straight filaments tend to be more luminous, more massive, 
having more massive central black holes and higher star for- 
mation rate. 

To see which property is most strongly correlated with 
Rl, we calculate the linear correlation coefficient between 
Rl and Pq as: 

_ {RlPg) 
^''^ ^{ARl){AP^y 

Table [2] lists the values of ^p^, for the four different galaxy 
properties. As can be seen, the central black hole mass is 
most strongly correlated with Rl - The star formation rate is 
also quite strongly correlated with Rl- This results indicate 
that the the coherent matter-accretion and gas-flow along 
the straight fllaments feed the central black holes and trigger 
the star formation efficiently. 



^ A^total — Ai'star + Afbulgo + A^coldGas + A^hotGas + A^ejected + 

Mcbh 
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Figure 3. Means of the central black hole mass (top-left), the stellar mass (top-right), the total mass (bottom- left), and the star 
formation rate (bottom-right) averaged over the constituent galaxies as a function of the host filament's linearity. In each panel the 
horizontal dotted line corresponds to the mean values averaged over all filaments. 



Table 2. The linear correlation coefficients. 
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0.89 


0.69 0.67 


0.85 



3.2 Linearity vs. Halo Mass 

Now that strong correlations between tlie linearity of void 
filaments and the physical properties of void galaxies are 
found, we would like to investigate the dependence of the 



mass of void halos on the linearity of their filament's linear- 
ity, using the catalogue of the filaments of void halos that 
was constructed in PL09. For each filament of void halos, we 
determine the maximum mass {Mmax) and the minimum 
mass [Mmin) among the masses of the constituent halos. 
Then, we calculate the means of Mmax and Mmin as a func- 
tion of Rl- 

Fig. |4] plots {Mmax) and {Mmin) as a function of Rl 
in the left and right panels, respectively. The errors in each 
panel is calculated as (AM^)^''^/v'A^. As can be seen, the 
value of {Mmax) decreases almost monotonically with Rl, 
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while the value of (Mmin) increases monotonically with Rl- 
It indicates that the differences among the masses of void ha- 
los are smaller when the parent filaments are more straight. 
To see this more clearly, we calculate the mean values of 
Mjnax — Mmin sls a function of Rl which is plotted in 
Fig. [5] It is now obvious that the difference decreases as 
Rl increases. In other words, the dark halos that constitute 
straight filaments tend to have similar masses. The efficient 
matter-accretion along straight filaments tend to synchro- 
nize the masses of void halos and regulate the shapes of 
potential wells of the void filaments to be flat. 



4 SUMMARY AND DISCUSSION 

We have found strong correlations of the stellar mass, cen- 
tral black hole mass and star formation rate of void galaxies 
with the linearity of their parent filaments. Among the three 
properties of void galaxies, the masses of the central black 
holes are mostly strongly correlated with the linearity of void 
filaments. The existence of this cross-correlations indicates 
the bridge effect of void filaments. Matter and gas from the 
surrounding large-scale structures accrete and fiow onto the 
void galaxies more efficiently along the more straight fila- 
ments. The dark halos constituting straight filaments are 
found to have similar masses, which implies that the poten- 
tial wells generated by straight void filaments have rather 
fiat-shaped bottoms. The efficient accretion of matter along 
straight filaments tend to reduce mass difference among the 
constituent void halos. 

It is interesting to note that our results may pro- 
vide an explanation on high specific star formation rate 
of void galaxies and high activity of voi d AGNs found in 
data from the Sloan Digital Sky Survey (|Constantin et al] 




0.2 0.4 0.6 0.8 1 

Figure 5. Means of the differences between the maximum and 
minimum masses of void halos as a function of the linearity of 
their parent's filaments, Rl- 



l2008l : lRoias et al.ll2005l '). The void filaments should be more 
straight on average than the wall filaments since the non- 
linear effect that can distort the filaments is less strong 
in void regions. Thus, matter-accretion and gas-fiow along 
more straight void filaments are likely to supply more effi- 
ciently fuel for the stars and AGNs of the void galaxies. 
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It is, however, worth noting that our results are sub- 
ject to the specific choice of the void-finding algorithm and 
the filament-finding algorithm. Unlike bound halos, neither 
voids nor filaments can be uniquely defined. A s shown in 
the current seminal work of lColberg et al.l (|2008| ). the voids 
identified by the HV02 algorithm tend to be relatively large 
in size compared with those found by different void-finders. 
Thus, if other algorithms were used, then it would have been 
difficult to find long filaments in voids. 

Besides, our results are based on numerical data from 
the Millennium Run simulations. To confirm the existence 
of the bridge effect of void filaments, an observational anal- 
ysis dealing with real data is required. Unfortunately, how- 
ever, it is quite unlikely to obtain any meaningful statistics 
about the bridge effect of void filaments from the c urrently 
available observational data. As first pointed out bv lPeeblei 
l|200ll ). the observed voids look apparently much emptier 
than the simulated vo ids, containing very few galaxies (e.g., 
iFoster fc NelsonI 120091 ). It is believed that this discrepancy 
between simulation and observation on voids indicates not 
the failure of the concordance cosmology but the limitation 
of current observational technique. Anyway it would be ex- 
tremely difficult to find as many void filaments as necessary 
for statistical analysis from the current available observa- 
tional data. Hence, we think that an observational test of 
our results is beyond the scope of this paper. It is concluded 
that our numerical finding of the bridge effect of void fila- 
ments may provide a new insight about how galaxies form 
and evolve in a ACDM cosmology. 
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